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We study experimentally and theoretically transition phenomena caused by the external forcing 
from Rayleigh-Benard convection with the large-scale circulation (LSC) to the limiting regime of 
unstably stratified turbulent flow without LSC whereby the temperature field behaves like a passive 
scalar. In the experiments we use the Rayleigh-Benard apparatus with an additional source of 
turbulence produced by two oscillating grids located nearby the side walls of the chamber. When 
the frequency of the grid oscillations is larger than 2 Hz, the large-scale circulation (LSC) in turbulent 
convection is destroyed, and the destruction of the LSC is accompanied by a strong change of the 
mean temperature distribution. However, in all regimes of the unstably stratified turbulent flow 
the ratio [(^.V.r)^ -|- (lyVyTf + ^Tf] / {6^) varies slightly (even in the range of parameters 
whereby the behaviour of the temperature field is different from that of the passive scalar). Here li 
are the integral scales of turbulence along x, y, z directions, T and 9 are the mean and fluctuating 
parts of the fluid temperature. At all frequencies of the grid oscillations we have detected the 
long-term nonlinear oscillations of the mean temperature. The theoretical predictions based on the 
budget equations for turbulent kinetic energy, turbulent temperature fluctuations and turbulent 
heat flux, are in agreement with the experimental results. 

PACS numbers: 47.27.te, 47.27.-i 



I. INTRODUCTION 

Various aspects of formation of large-scale circulations 
(LSC) and other types of coherent structures in a tur- 
bulent convection at large Rayleigh numbers have been 
discussed in a number of studies for the atmospheric flows 
[ll-Q , laboratory experiments in the Rayleigh-Benard ap- 
paratus 5-,30j and direct numerical simulations [3ll434l [. 
A mean-field theory [1, IH, HI] has been proposed for the 
explanation of formation of coherent structures in turbu- 
lent convection. It is known from laboratory experiments 
and numerical simulations 0, [11, [IE [U that the coher- 
ent structures in the Rayleigh-Benard turbulent convec- 
tion are not driven by the turbulent Reynolds stresses, 
associated with the tilting plumes at the upper and the 
lower horizontal walls. If the mean flow is established, 
the temperature of the fluid is larger at one side of the 
LSC and smaller at the other side of the LSC, and the 
mean flow is driven by the mean buoyancy force at the 
side walls of the chamber. 

The effect of coherent structures (mean wind) on global 
properties of turbulent convection is a subject of numer- 
ous discussions (see, e.g., reviews [37l - [39f ). It has been 
recently found in laboratory experiments in (40| . that 
in turbulent convection with the mean wind the bulk 
production of turbulence is due to the large-scale shear 
caused by the mean wind rather than by buoyancy. 
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Numerous experiments demonstrate that the mean 
temperature distribution in turbulent convection is 
strongly inhomogeneous and anisotropic due to the mean 
wind (see, e.g., [3, HI, [ll [11 [i^l ) , and the mean tem- 
perature gradient in the horizontal direction inside the 
LSC can be significantly larger than in the vertical di- 
rection, i.e., the vertical turbulent heat flux inside the 
LSC is very small. The dependencies of the dissipation 
rate (normalized by the molecular temperature diffusion) 
of the temperature fluctuations {0"^) as a function of the 
Rayleigh number and different statistical characteristics 
of this normalized dissipation rate have been studied ex- 
|)erimentally for turbulent Rayleigh-Benard convection in 

On the other hand, in a forced turbulent convection 
without a mean wind (i.e., in unstably stratified turbu- 
lent flows) the temperature field behaves like a passive 
scalar (see, e.g., |37l. l4ll - [43| |). It is not clear why in tur- 
bulent convection with LSC the vertical gradient of the 
mean temperature field is small inside LSC, and the verti- 
cal temperature gradient can even change the direction, 
i.e., it can be directed upwards. A detailed study of a 
transition from one regime with the mean wind to an- 
other regime without the mean wind due to additional 
forcing in unstably stratified turbulent flow, can elucidate 
the effects of the mean wind on turbulent convection. 

Different effects in turbulent convection with a forcing 
have been studied previously (see, e.g., fill, [d^. [46l [47|). 
In particular, an experimental study of a turbulent con- 
vection with forcing in the form of periodical pulses has 
been performed in j44| . In these experiments, in spite 
of the forcing, a large-scale circulatory flow was not de- 
stroyed, and the heat transfer in this system was en- 
hanced in comparison with constant and sinusoidally 
modulated energy inputs. It was found in {3] that the 



enhancement of the heat transfer depends on the syn- 
chronization of the kicking period of energy input with 
the intrinsic time scale of the turbulent flow. In other 
type of experiments [i^, H^l a relation between convec- 
tive heat flux and temperature gradient was determined 
in a vertical channel fllled with water, the average vertical 
mass flux being zero. It was shown in these experiments 
that when the LSC is suppressed by physically blocking 
its path, there exists a strong vertical temperature gradi- 
ent in the bulk of the Rayleigh-Benard apparatus and the 
temperature remains active, rather than passive field. 

The goal of this paper is to study experimentally and 
theoretically transition phenomena due to the external 
forcing from Rayleigh-Benard convection with LSC to 
the limiting regime of unstably stratified turbulent flow 
without LSC, whereby the temperature held behaves like 
a passive scalar. In the experimental study we use the 
Rayleigh-Benard apparatus with an additional source of 
turbulence produced by the two oscillating grids located 
nearby the side walls of the chamber. Additional forc- 
ing for turbulence allows us to observe evolution of the 
mean temperature and velocity fields during the transi- 
tion from turbulent convection with the LSC (for very 
small frequencies of the grid oscillations) to the limiting 
regime of unstably stratified flow without LSC (for very 
high frequencies of the grid oscillations). In the exper- 
iments we use Particle Image Velocimetry to determine 
the turbulent and mean velocity flelds, and a specially de- 
signed temperature probe with sensitive thermocouples is 
employed to measure the temperature field. 

Note that this experimental set-up with the turbulent 
convection and the external forcing has been previously 
used in studies of turbulent transport of particles (see, 
e.g., [isl- lisl - lsoj ). Comprehensive investigation of turbu- 
lence structures, mean velocity and temperature distribu- 
tions, velocity and temperature fiuctuations can elucidate 
a complicated physics related to particle clustering and 
formation of large-scale inhomogeneities in particle spa- 
tial distributions in unstably stratified turbulent flows. 
Notably, turbulent convection with the external forcing 
is used in chemical industrial devices to improve mixing 
(see, e.g., (FTl - tSBj ). Another example of turbulent con- 
vection with forcing is atmospheric wind superimposed 
on turbulent convective boundary layer. 

The paper is organized as follows. Section II describes 
the experimental set-up, instrumentation and the results 
of laboratory study of the transition phenomena in un- 
stably stratifled turbulent flow. The theoretical analysis 
and comparison with the experimental results are per- 
formed in Section III. Finally, conclusions are drawn in 
Section IV. In Appendix we present the detailed deriva- 
tions of the theoretical model for a passive scalar that 
used for the explanation of the observed effects. 
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FIG. 1. Experimental set-up: (1) -laser light sheet; (2) - heat 
exchangers; (3) - grid driver; (4) - digital CCD camera. 



II. EXPERIMENTAL STUDY 

Let us start with description of the experimental set- 
up and the results of laboratory study of the transition 
phenomena in unstably stratifled turbulent flow. 



A. Experimental set-up and instrumentation 

The experiments on transition phenomena in unstably 
stratifled turbulence have been conducted in rectangu- 
lar chamber with dimensions 26 x 58 x 26 cm^ in air 
flow with the Prandtl number Pr — 0.71. The side walls 
of the chambers are made of transparent Perspex with 
the thickness of 10 mm. In the experiments we use the 
Rayleigh-Benard apparatus with an additional source of 
turbulence produced by two oscillating grids. Pairs of 
vertically oriented grids with bars arranged in a square 
array (with a mesh size 5 cm) are attached to the right 
and left horizontal rods (see Fig. [1]). The grids are posi- 
tioned at a distance of two grid meshes from the chamber 
walls and are parallel to the side walls. Both grids are 
operated at the same amplitude of 61 mm, at a random 
phase and at the same frequency which is varied in the 
range from 0.5 Hz to 16.5 Hz. Here we use the following 
system of coordinates: z is the vertical axis, the y-axis is 
perpendicular to the grids and the xz-plane is parallel to 
the grids. 

A vertical mean temperature gradient in the turbulent 
air flow was formed by attaching two aluminium heat ex- 
changers to the bottom and top walls of the test section 
(a heated bottom and a cooled top wall of the chamber). 
In order to improve heat transfer in the boundary lay- 
ers at the bottom and top walls we used heat exchangers 
with rectangular flns 3 x 3 x 15 mm'^ (see Fig.[T]) which al- 
lowed us to form a mean temperature gradient up to 115 
K/m at a mean temperature of about 308 K when the 
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frequency of the grid oscillations / > 10 Hz. The thick- 
ness of the aluminium plates with the fins is 2.5 cm. The 
top plate is a bottom wall of the tank with cooling water. 
Temperature of water circulating through the tank and 
the chiller is kept constant within 0.1 K. Cold water is 
pumped into the cooling system through two inlets and 
flows out through two outlets located at the side wall 
of the cooling system. The bottom plate is attached to 
the electrical heater that provides constant and uniform 
heating. The voltage of a stable power supply applied to 
the heater varies up to 200 V. The power of the heater 
varies up to 300 W. 

The temperature field was measured with a tem- 
perature probe equipped with twelve E-thermocouples 
(with the diameter of 0.13 mm and the sensitivity of 
« 65 /iV/K) attached to a vertical rod with a diameter 4 
mm. The spacing between thermal couples along the rod 
was 22 mm. Each thermocouple was inserted into a 1 mm 
diameter and 45 mm long case. A tip of a thermocouple 
protruded at the length of 15 mm out of the case. The 
temperature was measured for 12 rod positions with 23 
mm intervals in the horizontal direction, i.e., at 144 lo- 
cations in a flow. The exact position of each thermocou- 
ple was measured using images captured with the optical 
system employed in PIV measurements. A sequence of 
600 temperature readings for every thermocouple at ev- 
ery rod position was recorded and processed using the 
developed software based on Lab View 7.0. Temperature 
maps were plotted using Matlab 7. 

In order to decrease mean flow generated by oscil- 
lating grids and to increase the range of homogeneous 
and isotropic turbulence we use partitions located in xz 
plane between grids and side walls of the chamber (of 
the width 85 mm). The aspect ratio of the chamber 
A = Hy/Hz — 1.52, where Hy is the size of the chamber 
along y-axis between partitions and is the hight of 
the chamber z-axis, respectively. 

The velocity fields were measured using a Stereoscopic 
Particle Image Velocimetry (PIV), see [54-56]. In the 
experiments we used LaVision Flow Master III system. 
A double-pulsed light sheet was provided by a Nd-YAG 
laser (Continuum Surelite 2 x 170 mJ). The light sheet 
optics includes spherical and cylindrical Galilei telescopes 
with tuneable divergence and adjustable focus length. 
We used two progressive-scan 12 bit digital CCD cam- 
eras (with pixel size 6.7 /im x6.7/im and 1280 x 1024 
pixels) with a dual-franie-technique for cross-correlation 
processing of captured images. A programmable Timing 
Unit (PC interface card) generated sequences of pulses 
to control the laser, camera and data acquisition rate. 
The software package LaVision DaVis 7 was applied to 
control all hardware components and for 32 bit image ac- 
quisition and visualization. This software package com- 
prises PIV software for calculating the flow fields using 
cross-correlation analysis. 

In order to obtain velocity maps in the central region 
of the flow in the cross-section parallel to the grids and 
perpendicular to a front view plane, we used one cam- 



era with a single-axis Scheimpflug adapter. The angle 
between the optical axis of the camera and the front 
view plane as well as the angle between the optical axis 
and the probed cross-section was approximately 45 de- 
grees. The perspective distortion was compensated us- 
ing Stereoscopic PIV system calibration kit whereby the 
correction was calculated for a recorded image of a cal- 
ibration plate. The corrections were introduced in the 
recorded in the probed cross-section images before their 
processing using a cross-correlation technique for deter- 
mining velocity fields. 

An incense smoke with sub-micron particles {pp/ p ^ 
lO'^), as a tracer was used for the PIV measurements. 
Smoke was produced by high temperature sublimation 
of solid incense grains. Analysis of smoke particles us- 
ing a microscope (Nikon, Epiphot with an amplification 
of 560) and a PM-300 portable laser particulate analyzer 
showed that these particles have an approximately spher- 
ical shape and that their mean diameter is of the order 
of 0.7/im. The probability density function of the par- 
ticle size measured with the PM300 particulate analyzer 
was independent of the location in the flow for incense 
particle size of 0.5 — 1 /im. The maximum tracer particle 
displacement in the experiment was of the order of 1 /4 
of the interrogation window. The average displacement 
of tracer particles was of the order of 2.5 pixels. The av- 
erage accuracy of the velocity measurements was of the 
order of 4% for the accuracy of the correlation peak de- 
tection in the inte rrog ation window of the order of 0.1 
pixel (see, e.g., [5J-[5a])- 

We determined the mean and the r.m.s. velocities, two- 
point correlation functions and an integral scale of tur- 
bulence from the measured velocity fields. Series of 520 
pairs of images acquired with a frequency of 1 Hz, were 
stored for calculating velocity maps and for ensemble and 
spatial averaging of turbulence characteristics. The cen- 
ter of the measurement region in yz-plane coincides with 
the center of the chamber. We measured velocity in a 
flow domain 256 x 256 mm^ with a spatial resolution of 
1024 x 1024 pixels. This corresponds to a spatial reso- 
lution 250 fiui I pixel. The velocity field in the probed 
region was analyzed with interrogation windows of 32 x 32 
or 16 X 16 pixels. In every interrogation window a veloc- 
ity vector was determined from which velocity maps com- 
prising 32 X 32 or 64 x 64 vectors were constructed. The 
mean and r.m.s. velocities for every point of a velocity 
map were calculated by averaging over 520 independent 
maps, and then they were averaged over the central flow 
region. 

The two-point correlation functions of the velocity field 
were determined for every point of the central part of the 
velocity map (with 16 x 16 vectors) by averaging over 
520 independent velocity maps, which yields 16 correla- 
tion functions in y and z directions. Then the two-point 
correlation function was obtained by averaging over the 
ensemble of these correlation functions. An integral scale 
of turbulence, £, was determined from the two-point cor- 
relation functions of the velocity field. In the experiments 
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we evaluated the variability between the first and the last 
20 velocity maps of the series of the measured velocity 
field. Since very small variability was found, these tests 
showed that 520 velocity maps contain enough data to 
obtain reliable statistical estimates. We found also that 
the measured mean velocity field is stationary. 

The characteristic turbulence time in the experiments 
Tz = 0.28 — 0.62 seconds, while the characteristic time pe- 
riod for the large-scale circulatory fiow 10 s) is by 1.5 
orders of magnitude larger than . These two character- 
istic times are much smaller than the time during which 
the velocity fields are measured 600 s). The size of 
the probed region did not affect our results. The temper- 
ature difference between the top and bottom plates, AT, 
in all experiments was 50 K (i.e., the global Rayleigh 
number based on molecular transport coefficients, was 
Ra 0.73 X 10^). Similar experimental set-up and data 
processing procedure were used previously in the experi- 
mental study of different aspects of turbulent convection 
[13,1101 9.nd in [isllisl - fsol for investigating a phenomenon 
of turbulent thermal diffusion [53, Si • 



B. Experimental results 

Let us start the analysis of the obtained experimental 
results from the dynamics of the mean velocity field. In 
Fig. [5] we show the mean flow patterns obtained in the 
experiments with unforced turbulent convection (Fig. [5^) 
and with forced turbulent convection at different frequen- 
cies of the grid oscillations (Figs. HJa anddJ:). The solid 
lines in Fig. [5] are stream lines in the cental yz-plane for 
X = 13 cm. The center of Figs. [5^, [Ud, Ht, which has 
coordinates y = 10 cm and z = 10 cm, coincides with 
the center of the chamber. When the frequency of the 
grid oscillations is larger than 2 Hz, the large-scale cir- 
culation (LSC) in turbulent convection in the yz-plane is 
not observed (see Figs.[2}D and[2j;). On the other hand, a 
complicated mean flow of the spiral form in the xz-plane 
with a non-zero mean vorticity in the y direction exists 
when the frequencies of the grid oscillations 2 < / < 7 
Hz. 

The destruction of the LSC is accompanied by a strong 
change of the mean temperature distribution. The distri- 
bution of the measured mean temperature field depends 
strongly on the frequency / of the grid oscillations (see 
Figs. [3][S]). In particular, for the very low frequency / 
the thermal structure inside the LSC in turbulent con- 
vection is inhomogeneous and anisotropic. The hot ther- 
mal plumes concentrate at one side of LSC, and cold 
plumes accumulate at the opposite side of LSC [H, H^j . 
In the central part of the flow the vertical mean temper- 
ature gradient VzT changes its sign depending on the 
frequency of the grid oscfllations (see Fig. [3]). 

At all frequencies of the grid oscillations we have ob- 
served the long-term nonlinear oscillations of the mean 
temperature with the period that is of the order of 10 
s. Note that the temperature was measured at 144 lo- 




FIG. 2. Mean flow patterns obtained in the experiments with 
unforced turbulent convection (upper panel) and with forced 
turbulent convection at the frequencies / = 2.2 Hz (middle 
panel) and / = 16.5 Hz (lower panel) of the grid oscillations 
for the unstably stratified turbulent flow at the temperature 
difference between the top and bottom walls AT = 50 K. 
Coordinates y and z are measured in mm. 
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FIG. 3. Vertical gradient of the mean temperature VzT 
(squares) and the amplitude of the nonlinear long-term os- 
cillations of this gradient (triangles) versus the frequency / of 
the grid oscillations for the unstably stratified turbulent flow. 
The mean temperature gradient is measured in K cm~^ and 
the frequency / is measured in Hz. 
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FIG. 4. Horizontal gradient of the mean temperature VyT 
(squares) and the amplitude of the nonlinear long-term oscil- 
lations of this gradient (triangles) versus the frequency / of 
the grid oscillations for the unstably stratified turbulent flow. 
The mean temperature gradient is measured in K cm~^ and 
the frequency / is measured in Hz. 



cations in a flow (the spacing between thermal couples 
was 22 mm). We perform the sliding averaging of the 
instantaneous temperature field over time of 3 s (the 
characteristic turbulence time in the experiments was 
0.3 - 0.6 s) in order to determine the mean temperature 
field T at 144 locations. Then we determine the long- 
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FIG. 5. Horizontal gradient of the mean temperature VxT 
(squares) and the amplitude of the nonlinear long-term oscil- 
lations of this gradient (triangles) versus the frequency / of 
the grid oscillations for the unstably stratified turbulent flow. 
The mean temperature gradient is measured in K cm~^ and 
the frequency / is measured in Hz. 



term variations of the vertical mean temperature gradi- 
ent S{ViT) = ViT—ViT due to the nonlinear oscillations 
of the mean temperature, where i — x,y, z. The bar here 
denotes additional time averaging. Note that the sep- 
aration distance of 22 mm between thermal couples is 
sufficient to measure the gradients of the mean temper- 
ature. Indeed, the integral scale of turbulence is about 2 
cm and the characteristic length scale of the mean tem- 
perature field is much larger than the integral scale of 
turbulence. Note also that for large frequencies of the 
grid oscillations the role of the thermal boundary layer 
in turbulent convection diminishes. It must be empha- 
sized that in all our experiments rectangular fins were 
attached at the bottom and top walls of the chamber in 
order to improve heat transfer in the boundary layers. 

In Fig. [6] we show time dependencies of the instan- 
taneous temperature T*°* = T + 9, the long-term vari- 
ations of mean temperature 6T — T ~ T and the long- 
term variations of the vertical mean temperature gradient 
^(VzT) = WzT — VzT due to the nonlinear oscillations 
of the mean temperature. In Figs. [3][5] we also show by 
triangles the amplitude, {[(5(ViT)]2}^/^, of the nonlinear 
long-term oscillations of the gradients of the mean tem- 
perature, where i = x,y, z. The dependencies shown in 
Figs, min] are obtained at the central region of the cham- 
ber with the size of 10 x 10 x 10 cm"^. 

In Fig. [7] we show the results of a Fourier analysis of the 
signal 6T = T — T. Inspection of Fig. [7] shows that there 
are two main maxima in the spectrum with the periods 
12.5 s and 20 s. Other smaller maxima in the spectrum 
are at the frequencies which are multiples of these main 
frequencies or their sums and differences. These are typ- 
ical features of nonlinear oscillations. The theory that 
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FIG. 6. Time dependencies of the instantaneous temperature 
rptot ^T + e, the variations of mean temperature ST ^ T -T 
and the variations of the vertical mean temperature gradient 
5(VzT) — V zT — VzT due to the long-term nonlinear oscil- 
lations of the mean temperature (with the period ~ 10 s). 
The bar denotes additional time averaging. These time de- 
pendencies are measured in the center of the chamber at the 
frequency / = 10 Hz of the grid oscillations for the unstably 
stratified turbulent flow. These temperature characteristics 
are measured in K. 
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FIG. 7. The normalized spectrum function E„{f) = 
l('5r)/lV I\iST)ffdf of the signal ST = T{t) - T, where 
in the Fourier space (<5T)/ = f ST exp[—ift] dt and / is the 
frequency of the nonlinear long-term oscillations of the mean 
temperature. 



explains the mechanism of these nonlinear oscillations of 
the mean temperature field, has not been developed yet. 
One may hypothesize that there are two possible mech- 
anisms for such oscillations. One mechanism could be 
related to the large-scale ToUmien-Schlichting waves in 
sheared turbulent flows (see [1^), which can cause the 
nonlinear oscillations of the mean temperature field. An- 
other mechanism of the nonlinear oscillations could be 
related to the generation of small-scale kinetic helicity 
due to large-scale shear flows in the system. The large- 
scale shear generates large-scale helicity, and since the 
total helicity is conserved, a non-zero small-scale helic- 
ity is produced. Our preliminary analysis shows that the 
small-scale helicity can cause the nonlinear oscillations of 
the mean temperature field. 

In order to avoid side effects of the grids we present the 
experimental results recorded in the central region of the 
chamber with the size of 10 x 10 x 10 cm'^. In Figs. 
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FIG. 8. Components u™" = \/(it|) (triangles) and u™" = 
(ti|) (squares) of the r.m.s. velocity versus the frequency 
/ of the grid oscillations for the unstably stratified turbulent 
flow. The turbulent velocity is measured in cm s~^ and the 
frequency / is measured in Hz. 



we show the frequency dependencies of the following mea- 
sured turbulence parameters: the r.m.s. velocity fluctu- 
ations ^ (uy) and ^/(ufy, and the corresponding integral 
scales of turbulence along horizontal y and vertical z di- 
rections {£y and £z), where u are the fluctuations of the 
fluid velocity. The mean, r.m.s. velocities and the two- 
point correlation functions of the velocity fleld have been 
calculated by averaging over 520 independent velocity 
maps, and then they have been averaged over the central 



flow region. The integral scales of turbulence, £y and iz, 
have been determined from the normalized two-point lon- 
gitudinal correlation functions of the velocity fleld, e.g., 
Pyiy) = {uy{ro)uy{ro+yey))/{ul{ro)) [and similarly for 
Fz {z) after replacement in the above formula y hy z], us- 
ing the following expression: £y — {J^ Fy{y)dy)s [and 
similarly for £z after replacement in the above formula y 
by z], where L = 10 cm is the linear size of the probed 
flow region, By is the unit vector in the y direction, and 
(...)s is the additional averaging over the yz plane. 

Inspection of Fig. [3] shows that the turbulent length 
scales are weakly dependent on the frequency / of the 
grid oscillations in the unstably stratified turbulent flow 
for / > 3 Hz, while the turbulent velocities and turbu- 
lent time scales vary strongly with the frequency /. Note 
that Ty f and / tend to a constant 3) for higher fre- 
quencies (/ > 5 Hz) of the forcing, where Ty — £y/ ^ (Uy) 

and Tz = £zl \/ (Uz) are the characteristic turbulent times 
in horizontal and vertical directions. Note also that the 
ratio of the characteristic values of the mean velocity to 
the turbulent velocity in the central part of the chamber 
varied from 1 at / = to 0.5 when the frequency / of 
the grid oscillations / > 2 Hz. In particular, the char- 
acteristic values of the mean velocity at the frequencies 
of the grid oscillations / = 0; 2.2; 16 Hz are as following: 
7.5 cm/s; 3 cm/s and 6 cm/s, respectively. 

In Fig. [TU] we show the frequency dependence of the 
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FIG. 9. Horizontal and vertical integral scales of turbulence 
£y (triangles) and £z (squares) versus the frequency / of the 
grid oscillations for the unstably stratified turbulent flow. The 
turbulent length scales are measured in mm and the frequency 
/ is measured in Hz. 



ratio u*/u™^, where u™* is the r.m.s. of the vertical 
component of velocity fluctuations in the unstably strat- 
ified turbulent flow and w* is the r.m.s. of the vertical 
component of velocity fluctuations in the isothermal tur- 
bulence. In Fig. [TU] along with the ratios of the r.m.s. of 
the vertical component of velocity fluctuations, it*/u™*, 
in isothermal and unstably stratified turbulent flows, we 
also plot the ratio ul/iiz, where Uz is the vertical com- 
ponent of the effective turbulent velocity, 



(1) 



P = g/Ti, is the buoyancy parameter, T* is a reference 
value of the mean absolute temperature and g is the ac- 
celeration of gravity. This definition of the effective ve- 
locity is derived from the budget equation (j8]) for the 
turbulent kinetic energy in Sect. III. This effective veloc- 
ity takes into account the production of the turbulence 
by buoyancy. Inspection of Fig. [10] shows that the val- 
ues of these ratios are very close. The latter implies that 
the measured turbulent velocity in the unstably stratified 
turbulent flow, it™*, is of the order of Uz- 

The dependence of the r.m.s. of the temperature fluc- 
tuations (0^) versus the frequency / of the grid oscilla- 
tions is shown in Fig. 1111 This dependence is caused by 
the nontrivial behaviour of the mean temperature gradi- 
ents along X, y, z axes (see Figs. |3]l5]), where Q are fluctu- 
ations of fluid temperature. We have also analyzed the 
experimental results obtained at the smaller temperature 
difference, AT = 25 K, between the top and bottom walls 
(i.e., smaller Ra), and found only minor changes in the 
final results. In the next section we analyze the obtained 
experimental results. 
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FIG. 10. Ratios ul/u™" (squares) and ul/uz (circles) ver- 
sus the frequency / of the grid oscillations. Here is 
the r.m.s. of the vertical component of velocity fluctuations 
in the unstably stratified turbulent flow, u* is the r.m.s. of 
the vertical component of velocity fluctuations in the isother- 
mal turbulence, Uz is the vertical component of the effective 
turbulent velocity, = [{Kf + 44 P that takes 

into account the production of the turbulence by buoyancy. 
The velocity is measured in cm s~^ and the frequency / is 
measured in Hz. 
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FIG. 11. The r.m.s. of temperature fluctuations 9rma = 
(9^) versus the frequency / of the grid oscillations for the 
unstably stratified turbulent flow. 



III. THEORETICAL ANALYSIS AND 
COMPARISON WITH EXPERIMENTAL 
RESULTS 

Let us first analyze the frequency dependence of the 
temperature fluctuations. To this end we consider the 
budget equation for (6^): 



Dt 



-div*e = -2(F-V)T-2ee 



(2) 



(see, e.g., [60|-l63|), where D/Dt = d/dt + U-V, u are 
the fluctuations of the fluid velocity, U is the mean ve- 
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locity that describes the coherent structures, = (u0^) 
is the third-order moment that determines the flux of 
(f?^), Fi — {ui9) is the turbulent heat flux, 9 are the 
temperature fluctuations, T is the mean temperature, 
ee = I?((V6')^) is the dissipation rate of (6*^), where D is 
the coefficient of molecular temperature diffusion. 

Let us use the Kolmogorov-Obukhov hypothesis which 
allows to estimate the dissipation rate eg of {0'^), as 
[d^)/To (see, e.g., [H, [el]), where tq = £/uo is 



the spectral r approximation that is valid for large Peclet 
numbers. 

Now let us discuss the generalization of Eq. (jlj. To 
this end we use the budget equation for the turbulent 
heat flux Fi — (ujO): 



DF\ 
Dt 



the characteristic turbulent time and Uq is the charac 
teristic turbulent velocity at the scale I. Indeed, eo 



(F-V)[/, - 



' Vip) - {uiUj) Vj T 
, (5) 



_ where /3i = 13 Cz, p are the pressure fluctuations, p is 



D{{VeY) = D{0^) j^'eEg{k)dk « (^?2)/ro, where ^^^^^ density, U is the mean velocity, is the 



Eg{k) — (q — l)kQ^ (fc/fco)~' is the spectrum function 
of the temperature fluctuations, q is the exponent of the 
spectrum of the temperature fluctuations (e.g., q = 5/3 
for the Corrsin-Obukhov spectrum), ko = , kb = i'^^, 
4 = £/Pe^/(^"«' and Pe = uo£/D is the Peclet number. 
The latter estimate implies that the main contribution 
to the dissipation rate eg arises from very small molec- 
ular temperature diffusion scales £},. The Kolmogorov- 
Obukhov hypothesis has been widely used in atmospheric 
turbulence in a number of studies including investigations 
of atmospheric stably stratified and convective boundary 
layers. In a steady-state Eq. 1^ yields: 



') « -2to (f- v)r 



(3) 



In deriving Eq. ([3]) we take into account that the com- 
ponents of the turbulent heat flux and turbulent temper- 
ature diffusion coefficients in an anisotropic turbulence 
are given by = V,T and = C^i, /(^, 

and similarly for other components Fy , F^ of the turbu- 
lent heat flux and the turbulent temperature diffusion 
coefficients Dy , _D J with the replacement in the above 
formulas x by ?/ or by z. We also neglected a small 
term div ^g in Eq. ([2]) for nearly homogeneous turbu- 
lent convection. When horizontal gradients of the mean 
temperature are much smaller than the vertical gradi- 
ents, |V;j,yT| -C |V^r|, Eq. © yields the following non- 
dimensional ratio 



const 



(4) 



Kraichnan (1968) derived Eq. Q for a passive scalar 
regime using the Lagrangian-history direct-interaction 
approximation for small Peclet numbers (66j . Equa- 
tion (HI) is in agreement with the corresponding equation 
derived by means of the path integral approach in [67 1 
for the passive scalar advected by the delta-correlated in 
time Gaussian smooth velocity field. Equation (j4]) can 
be also derived using spectrum of temperature fiuctua- 
tions found by dimensional ar gum ents in and by the 
renormalization procedure in 16911 . Appendix we also 
derive Eqs. Q and (jH) using other procedures, namely a 
quasi-linear approach (or a second order correlation ap- 
proximation) that is valid for small Peclet numbers and 



( PI 

dissipation rate of the turbulent heat fiux and = 
{uiUjO) -f 5ijp~^{9p)/2 is the term that includes the 
third-order moments. According to the estimate made 
in 62], /3, (6l2) -Hp-i {OV.p) « 2/3, (6*2). In a steady-state 
case Eqs. (O and ([5]) yield: 



1 



ICgCp 



const 



(6) 



where 



[4 V,T]2 = [(4V,T)2 ^ (^^Vyr)2 + (4V,T)2] 

y^[l +^CgCFpTl{V.T)Y\ (7) 



and Cp is an empirical constant. Here we take into 
account that the dissipation rate of the turbulent heat 

fiux is e^' = FjJCFTQ and k, Ty ^ = tq. 
In Eqs. 



and ([7]) we also have neglected the terms 
0{t I Lip\l^ I {L\.Lu)\, where Lt and Lu are the char- 
acteristic spatial scales of the mean temperature and ve- 
locity field variations. Equation ([6]) is valid in a general 
case including the passive scalar regime. In the case of 
the passive scalar regime whereby there is only the verti- 
cal gradient of the mean temperature, Eq. (jB)) coincides 
with Eq. dH). 

In Fig . [T2] we plot the non-dimensional ratio 
^* V*T/a/ (IP) versus the frequency / of the grid oscilla- 
tions for the unstably stratified turbulent fiow obtained 
in our experiments, where we assumed that — iy and 
CgCp = 1/12. Inspection of Fig. [12] shows that this 
non-dimensional ratio varies slightly even in the range 
of parameters whereby the behaviour of the temperature 
field is different from that of the passive scalar. Small 
deviations of the experimental results from the theoreti- 
cal predictions [see Eq. ©] may be caused by a non-zero 
term, div ^g and deviations from the steady-state due to 
the long-term oscillations of the mean temperature. 

Now let us derive Eq. ([T]) for the measured turbulent 
velocity in the unstably stratified turbulent fiow. We 
use the budget equation for the turbulent kinetic energy 
Ek = (u2)/2: 



DEk 
Dt 



- div = ~{ui uj) VjUi + l3Fz + (u-ff) - Sk, 



(8) 

where = p~^{up) + (1/2) (uu^) is the term that in- 
cludes the third-order moments, (u-fy) is the production 
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FIG. 12. The non-dimensional ratio V ,T /Orms versus the 
frequency / of the grid oscillations for the unstably stratified 
turbulent flow, where 9rms = \/ {0'^)- 




FIG. 13. The non-dimensional ratio \Ep\/Ek of the potential 
to kinetic energies versus the frequency / of the grid oscilla- 
tions for the unstably stratified turbulent fiow. The horizontal 
line corresponds to \Ep\/Ek ~ 0.1. 



of turbulence caused by the grid oscillations and Ek is the 
dissipation rate of the turbulent kinetic energy. In the 
steady-state Eq. dH) yields (u^) = 2ro [i^^ 5^ + (u-f/) -f- 
P Fg:], where is the turbulent viscosity. Let us intro- 
duce the characteristic velocity for isothermal turbulence 
(m*)^ = 2to [v^ Si + (u-f/)], where S"* is the large-scale 
shear that can appear in isothermal oscillating grid tur- 
bulence. Note that for large frequencies of the grid os- 
cillations the large-scale circulations caused by turbulent 
convection disappear and S* — >■ S"* . Since the vertical tur- 
bulent heat flux = —Dt^zT and for large frequencies 
of the grid oscillations |V^T| ~ [see Eq. (g])], 

we arrive at the following equation: 

u^r ^[{<f + ^^zPVW)Y'\ (9) 

where we take into account that the r.m.s. of the hori- 
zontal component of velocity fluctuations in the unsta- 
bly stratified turbulent flow (u^) ~ u*^ and u*^ = 
[(u*Y — (u^)^]^/^. This is the reason that the frequency 
dependencies of the ratios u^ju™^ and u%/uz obtained 
in our experiments (see Fig. [TU|) . are close. 

Using Eq. ([2]) we obtain the evolutionary equation for 
the turbulent potential energy Ep = {l3^ /2N^) {6^): 

DE 

+ div *p = Pp - PFz - ep, (10) 

(see, e.g., [fii, [g!]), where A^^ ^ /3V^T, *p = 
lp^/2N^)^e, Pp = -(/3VA^2)(F,,-V)T is the source 
(or sink) of the turbulent potential energy caused by the 
horizontal turbulent heat flux F/j = (uhO), Uh is the 
horizontal component of the velocity fluctuations and 
Ep = {(3^/2N'^)se. The turbulent potential energy is 
analogous to the available potential energy introduced by 
Lorenz [t^]: the available potential energy and the tur- 
bulent potential energy are proportional to the squared 
temperature fluctuations. The principal difference be- 
tween these two quantities is that the available potential 
energy is an integral property of the entire flow-domain 



(e.g., of the atmosphere as a whole), whereas the tur- 
bulent potential energy is a local quantity, i.e., it is de- 
termined in each point of turbulent flow. The buoyancy 
term, f3Fz, appears in Eqs. ([8]) and (fTO| with opposite 
signs and describes the energy exchange between the tur- 
bulent kinetic energy and the turbulent potential energy. 
These two terms cancel in the budget equation for the 
total turbulent energy, E — E^ + Ep-. 

DE 

— + div * = - {u, uj) VjU^ + (u-f/) - e, (11) 

where $ = + ^*p and e = £k + ep. The concept of the 
total turbulent energy is very useful in analysis of strat- 
ified turbulent flows. In particular, it allows to explain 
the physical mechanism for existence of turbulence for 
arbitrary values of the Richardson number, and abolish 
the paradigm of the critical Richardson number in the 
stably stratified atmospheric turbulence (see [6^ [63|). 

For very small horizontal mean temperature gradients, 
the source term Pp is much smaller than other terms 
in the right hand side of Eq. ([TT]). In Figs. [13] and [H 
we show the non-dimensional ratios \Ep\/ E^ and |Pp|/£fe 
versus the frequency / of the grid oscillations for the 
unstably stratified turbulent fiow. Inspection of Figs. [T3] 
and[T4]shows that the temperature fluctuations affect the 
turbulent kinetic energy only in the vicinities of / = 
and / = 5.2 Hz. Note that in the vicinity of / = 5.2 Hz, 
the sign of V zT changes from positive to negative (see 
Fig. [3]). 

In view of the above analysis, we consider the passive- 
like scalar behaviour of the temperature fleld in the kine- 
matic sense. In particular, if the temperature fluctua- 
tions (6*^) do not affect the turbulent kinetic energy, the 
temperature fleld can be considered as a passive scalar. 
In this case the evolution of the temperature fleld in 
a given turbulent velocity field is a kinematic problem, 
whereby there is no dynamic coupling between the tem- 
perature fiuctuations, (0^), and the turbulent kinetic en- 
ergy in spite of a non-zero turbulent heat flux, i.e., non- 
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FIG. 14. The non-dimensional ratio |Ppj/efc of the source (or 
sink) of the turbulent potential energy caused by the horizon- 
tal turbulent heat flux, to the dissipation rate of the turbulent 
kinetic energy versus the frequency / of the grid oscillations 
for the unstably stratified turbulent flow. The horizontal line 
corresponds to \Pp\/Ek ~ 1. 



zero correlations of the velocity and temperature fluctu- 
ations. When the effect of the temperature fluctuations 
on the turbulent kinetic energy cannot be neglected, the 
temperature is considered as an active field. This un- 
derstanding of the passive or active behavior of the tem- 
perature field is different from that based on the scal- 
in^ behaviour of the temperature structure function (see 

Note that the turbulent convection studied in our ex- 
periments is non-Boussinesq, i.e., div u ^ 0. In the 
theoretical analysis we use the anelastic approximation, 
div (p u) = 0. It should be noted that the deviation from 
the Boussinesq approximation in our experiments is less 
than 10 % (i.e., spatial variations of the fluid density are 
5p/p <{).!). 



IV. CONCLUSIONS 

In the present study we have investigated transition 
phenomena caused by the external forcing from Rayleigh- 
Benard convection with LSC to the limiting regime of 
unstably stratified turbulent flow without LSC, whereby 
the temperature fleld behaves like a passive scalar. The 
external forcing of turbulence is produced by two oscil- 
lating grids. When the frequency of the grid oscillations 
is larger than 2 Hz, the large-scale circulation (LSC) in 
the j/z-plane in turbulent convection is destroyed. On the 
other hand, a complicated mean flow of the spiral form in 
the xz-plane with a non-zero mean vorticity in the y di- 
rection exists when the frequencies of the grid oscillations 
2 < / < 7 Hz. The destruction of the LSC is accompa- 
nied by a strong change of the mean temperature distri- 
bution. In the central part of the flow the vertical mean 
temperature gradient changes its direction depending on 
the frequency of the grid oscillations. We have shown 
that for all regimes of the unstably stratifled turbulent 



flow the ratio [{i^^V^Tf + {(■y'^yTf + {l^^V ^Tf]/ {6^) 
varies slightly even in the range of parameters whereby 
the behaviour of the temperature field is different from 
that of the passive scalar. The experimental results ob- 
tained in this study are in an agreement with the the- 
oretical predictions based on the budget equations for 
turbulent kinetic energy, turbulent potential energy and 
turbulent heat flux. 

We have also found that in the unstably stratifled tur- 
bulent flow the vertical component of the turbulent ve- 
locity fleld is determined by Eq. ([9|). The latter equation 
is derived from the budget equation for the turbulent ki- 
netic energy. This implies that the buoyancy heat flux 
contributes to the production of the turbulent kinetic en- 
ergy in turbulent convection in the absence of the LSC. 
However, the turbulent kinetic energy in large Rayleigh 
number turbulent convection with coherent structures 
(LSC) is produced by shear, rather than by buoyancy 
(see (Si). 

The long-term nonlinear oscillations of the mean tem- 
perature have been detected for all regimes, from the tur- 
bulent convection with LSC to the passive scalar regime, 
at all frequencies of the grid oscillations. 



Appendix A: Temperature fluctuations 

Evolution of the instantaneous temperature fleld 
rtot(i,r) in a turbulent flow is determined by the fol- 
lowing equation: 



dt 



(v- V) Ttot + (7 - 1) (V-v) Ttot = D V'Tu 



(Al) 



where rtot(^, r) = r-l-6', D is the coefficient of molecular 
temperature diffusion, 7 = Cp/cv is the ratio of specific 
heats, V is the velocity field that satisfies to continu- 
ity equation in anelastic approximation for a low-Mach- 
number flow: 

V-(pv) = 0. (A2) 

Combining Eq. (|Aip and (jA2p we obtain the following 
equation: 



dt 



(v- V) Ttot = D V Vtc 



(A3) 



where v = 7V and we take into account that for a low- 
Mach-number flow without imposed external pressure 
gradient Vp/p « -VTtot/Ttot. 

Averaging Eq. (jA3p over ensemble of turbulent velocity 
fleld we obtain the equation for the evolution of the mean 
temperature fleld T(t, r): 



dT 
'dt 



(U-V)r+ V-F = Z? V^T 



(A4) 



where U = 7U, U is the mean fluid velocity, F (uS) is 
the heat flux, u = 7U and u are the velocity fluctuations. 
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Temperature fluctuations for small Peclet 
numbers 



In order to study temperature fluctuations for small 
Peclet numbers we use a quasi-linear approach or a 
second order correlation approximation (see, e.g., (zll)- 
Subtracting Eq. (jA4|) from Eq. (|A3|) yields equation for 
the temperature fluctuations: 



at 



(A5) 



where / — — (u-V)r is the source term and Q — 
V- [vlO — {uO)] + 9{u-'V)p/ p is the nonlinear term. In 
the expression for Q we neglected the term {9u)-V p/ p 
which is quadratic in large-scale spatial derivative. Let 
us neglect the nonlinear term but keep the molecular dif- 
fusion term in Eq. (IA5|) . For this reason this approach 
is called a quasi-linear or perturbation approach. This 
approximation for a given velocity field is valid only for 
small Peclet numbers (Pe <C 1), where Pe — UQljD. Let 
us rewrite Eq. (|A5[) in a Fourier space. Then the solution 
of Eq. (|A5|) is given by Q{ui, k) — G_D(aj, k)/(a;, k), where 
Gi3(cj,k) = (Dfc^ 

Let us apply a standard two-scale approach, whereby 
the non-instantaneous two-point second-order correlation 
function is written as follows: 



(Mi(ti,x) e'(i2,y)) y (wi(wi,ki)0(a;2,k2))exp[i(ki-x 
-|-k2-y) + 1(^1^1 + ^2^2)] duJ2 d\<ii d\<i2 
= j k) cxp[ik-r + ioj f ] dk , (A6) 

where we use large scale variables: R = (x -I- y)/2, K = 
ki -I- k2, t — {ti -\- ^2)72, VI = uji + W2, and small scale 
variables: r = x — y, k = (ki — k2)/2, f = ti — t2, 

CJ = (wi - W2)/2, 



Fi(w,k) = J (ui(wi,ki)6'(w2,k2)}exp[jr2t 

+iK-R]dndK, (A7) 

and uji =LLi + n/2, uj2 = -w + f2/2, ki = k + K/2, k2 = 
— k + K/2 (see, e.g., [z3)- We assume here that there 
exists a separation of scales, i.e., the maximum scale of 
random motions £ is much smaller than the characteristic 
scales of inhomogeneities of the mean temperature and 
fluid density. 

The turbulent heat flux Fi and the function Eg — 
(6'^)/2 are given by the following relations: 

Fi^ j {u.i{uj, k) I{-uj, -k)) du dk , (A8) 



x(v,r) (v,T). 



(A9) 



Hereafter we use the simple model for the second mo- 
ments of a random velocity field, (ui(a;, k) Mj(— w, — k)). 



in a Fourier space: 
(ui(w,k)wj(-a;,-k)) = 



l'<E{k) 



87rfc2 



fc2 



(AlO) 



where 5{uj) is the Dirac's delta-function, 5ij is the Kro- 
necker tensor, Pij{k) = Sij — kikj/k'^, the energy spec- 
trum function is E{k) = k^^ {q — 1) {k/ko)~'^, the expo- 
nent 1 < q < 3, the wave number fco = the length £ is 
the maximum scale of random motions and uq is the char- 
acteristic velocity in the maximum scale of random mo- 
tions. Note that contributions due to the density stratifi- 
cation in the random velocity field, {ui{uj, k) Uj{—uj, — k)) 
to the turbulent heat flux Fi and the temperature fluc- 
tuations (02) vanish. 

After integration in w-space and in k-space, Eqs. (|A8p 
and (|A9|1 yield formulae for the turbulent heat flux Fi and 
the temperature fluctuations {6'^) in an isotropic back- 
ground turbulence: 



(u, 9) = -DtV^T , Dt 
(6)2) = -C,f (VT)2, 



CduqI, 



(All) 
(A12) 



G.^^^^Pe, G, = i|^Pe^ 



3(g + i) 



3(g + 3) 



where Dt is the turbulent temperature diffusion coef- 
ficient. Equations (|Alip and (|A12[) also can be ob- 
tained using spectrum of temperature fluctuations for 
small Peclet numbers found by means of dimensional ar- 
guments in [73I and by the Lagrangian-history direct- 
interaction approximation applied in [66j . 



Temperature fluctuations for large Peclet 
numbers 



In this subsection we derive formulae for the turbulent 
heat flux Fi and the temperature fluctuation function 
Eg using the t approach that is valid for large Peclet 
and Reynolds numbers. Using Eq. (jA5[) written in a 
Fourier space we derive equation for the instantaneous 
two-point second-order correlation functions Fi{t,'k) ~ 
{u^{t, k) e{t, -k)) and Eg{t, k) = (1/2) (0(i, k) e{t, -k)): 



f]F 



MEj^''"{k), 



MF^^'"\k), (A13) 
(A14) 



where MF- 
72 



(///) 



(k) 



^[{u,Q) + {{dujdt)0) - 

D{u,Ve)]^ and ME^/"\k) = -(1/2) [(6'g) - 
D{d'V^9)]i^ are the third-order moment terms appearing 
due to the nonlinear terms which include also molecular 
diffusion term, and 

leit, k) = i [(/(t, k) 0(t, -k)) + {e{t, k) /(t, -k))] . 

The equation for the second moment includes the flrst- 
order spatial differential operators applied to the 
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third-order moments F^^^^K A problem arises how to 
close the system, i.e., how to express the third-order 
terms MF^^^^^ through the lower moments F'^^^ (see, 
e.g., [13, [H, [zl|)- We use the spectral r approxima- 
tion which postulates that the deviations of the third- 
moment terms, A4F'^^^^^(k), from the contributions to 
these terms afforded by the background turbulence, 
_A4F'^^^'''^(k), can be expressed through the similar de- 
viations of the second moments, F(^^)(k) — F(^^'°)(k): 

MF<'"'\k) -MF'-^^^'°\k) 
1 



F(")(k) -F(^^'°)(k) , (A15) 



(see, e.g., [7J-|76|), where Tr(fc) is the scale-dependent 
relaxation time, which can be identified with the corre- 
lation time r(fc) of the turbulent velocity field for large 
Reynolds and Peclet numbers. The functions with the 
superscript (0) correspond to the background turbulence 
with a zero gradient of the mean temperature. Validation 
of the T approximation for different situations has been 
performed in numerous numerical simulations and ana- 
lytical studies (see, e.g., review [73]; and also discussion 
in [7§\, Sec. 6). 

Note that the contributions of the terms with the 
superscript (0) vanish because when the gradient of 
the mean temperature is zero, the turbulent heat flux 
and the temperature fluctuations vanish. Consequently, 
Eq. (IMSl) for MF^"^\k) reduces to MFy"\k) = 

-F,(k)/T(fc) and for ME^/"\k) = -Ee{k)/T{k). We 
also assume that the characteristic time of variation of 
the second moments F'j(k) and Eg{k) are substantially 
larger than the correlation time T(fc) for all turbulence 
scales. Therefore, in a steady-state Eqs. (jA13l) and (IA14[) 
yield the following formulae for the turbulent heat flux 
Fj and the function Eg: 



Eg 



T{k){u,{t,k)I{t,-k))dk, (A16) 
T{k)Ig{t,k)dk. (A17) 



Now we use the following simple model for the second 
moments of turbulent velocity field, (ui(k) Uj{~k)) in k 
space: 



{ui{k) Uj{-k)) 



l^nlEjk) 
87rfc2 



fc2 



(A18) 



Contributions due to the fluid density stratification in 
the random velocity field, {ui{k) Uj{—k)) to the turbulent 
heat flux Fi and the temperature fluctuations (6'^) vanish. 
After integration in k-space Eqs. (jA16p and (jA17|) we 
arrive at equations for the turbulent heat flux, Fi, and the 
temperature fluctuations. Eg in an isotropic background 
turbulence: 

(u, 9} = ~DtS7,T , Dr^^Couoi, (A19) 

{0^) = -{8/9)fiVTf. (A20) 
where C/j = 7^/3. In the derivation we used the follow- 
ing expression for the turbulent correlation time, T(k) = 
2 To {k/koY~'', where tq = ^/mq is the characteristic tur- 
bulent time. Therefore, the formulae for the turbulent 
heat flux Fi and the function Eg are similar for small and 
large Peclet numbers, while the coefficients are different 
in these two cases. Notably, the theoretical predictions 
for large Peclet number (which correspond to the labora- 
tory experiments) are in agrement with experimental re- 
sults and with the estimate based on the balance equation 
for Eg [see Eq. ©]. Equations (|XT9|) and (1X20]) are in 
agreement with the corresponding equations derived by 
means of the path integral approach in [67| for the passive 
scalar advected by the delta-correlated in time Gaussian 
smooth velocity field. Equations (IA19P and (|A20p also 
can be obtained using spectrum of temperature fluctua- 
tions found by means of dimensional arguments in [6^ 
and by the renormalization procedure for large Peclet 
numbers performed in [69^. 

Note that in our experiments with large frequencies of 
the grid oscillations, i.e., in the limiting regime of the un- 
stably stratified turbulent flow without the LSC, whereby 
the temperature field behaves like a passive scalar, the 
constant Cd is nearly independent of the molecular heat 
transport coefficients for large Peclet numbers (which 
correspond to the conditions of our experiments). The 
constant Cd depends on the molecular heat transport 
coefficients only for small Peclet numbers. 
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